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ABSTRACT 

We explore the mass distribution of material associated with galaxies from the ob- 
servation of gravitational weak lensing for the galaxy mass correlation function with 
the aid of iV-body simulations of dark matter. The latter is employed to unfold various 
contributions that contribute to the integrated line of sight mass density. We conclude 
that galaxies have no definite edges of the matter distribution, extending to the middle 
to neighbouring galaxies with the density profile roughly r~ 2A beyond the virial radius. 
The mass distributed beyond the virial radius (gravitationally bound radius) explains 
the gap seen in the mass density estimates, the global value fi m ~ 0.27 and typically 
^gai ~ 0.15 from the luminosity density multiplied by the mass to light ratio. We sug- 
gest to use a physical method of gravitational lensing to characterise galaxy samples 
rather than characterise them with photometric means. 

Subject headings: 



1. Introduction 

Where and how is matter distributed in the Universe is a long-standing question. The widely 
accepted view is that matter of the universe is preponderantly borne by galaxies. This leads to 
an estimate of the mass density of the Universe (Ostriker et al. 1974). With modern parameters, 
the luminosity density L = 2.2 ± 0.3 x 10 8 /iL o Mpc" 3 (Blanton et al. 2001;2003) multiplied by 
the average mass to light ratio of gal axies, (M/L) ~ (17 ± 50)/iMo/L r Q, from the gravitational 



lensing shear around 2 x 10 4 galaxies ( McKay et al.ll200ll ) leads to the mass density of the Universe 
p m = £ x (M/L) as, 

n m = 0.13 ±0.05 (1) 

in units of the critical mass density. The mass of M/L we adopted here is supposed to represent 
the virial mass, i.e., the mass of particles gravitationally bound in galaxies. This M/L is consistent 
among a number of estimates for individual galaxies (e.g., Bahcall &; Fan 1998). 
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On the other hand, we now have a convincing estimate of the global matter density from the 



cosmic microwave background radiation anisotropy (jKomatsu et al.ll201ll ). which converges to 



a 



0.27 ±0.03. 



(2) 



The difference between Eqs. ([I]) and ([2]) raises the problem as to where is half the matter missing, 
indicated by the mismatch of the two numbers. 

Within the bound radius of galaxies, the mass distribution is inferred to obey the profile close 
to the one advocated by Navarro, Frenk & White (1997, hereafter NFW) in the CDM dominated 
universe. Simulations of clustering of dark matter conveniently point towards this distribution 
for mass assemblies. Th is is also supported by gravitational lensing analysis for bright galaxies 



(jMandelbaum et al.ll2006l ). Little is studied, however, concerning the mass distribution beyond the 
virial radius. 

Observational advancement with large galaxy samples lends gravitational weak lensing to a 
powerful tool to explore the average surface mass density of galaxies. The measure is the surface 
density S(i?) as a function of the projected radius from the centre of the galaxy R, 



pir 



(3) 



where p is the density at r and x is the line of sight distance, both measured from the centre of the 
galaxy. This quantity can be explored along the line-of-sight towards a light source by measuring 
its tangential shear 7 t or magnification p of the image of background sources, as 



fM(R) 



£(< R) - E(R) &L(R) 



1 + 2k(R) = 1 + 2 



E(R) 



(4) 
(5) 



where £(< R) is the average of £ within R, k is the convergence, and S cr = (c 2 / 'AttG)(D ' s / 'D ^ S D '^) 
is the critical surface density. Q With the modern large samples this is applicable even to the 
outskirt of galaxies significantly beyond the virial radius, say to ~ 10/i _1 Mpc using either lensing 
shear (Sheldon et al. 2004; Mandelbaum et al. 2006; Reyes et al. 2010) or statistical magnification 
(Menard et al. 2010). While thus obtained X(-R) includes the surface mass density from neighboring 
galaxies along the line of sight, this quantity gives a cogent information as to the surface mass 
density profile of the galaxy away from its central region. 



R- 
R 



Menard et al.l (|2010l . hereinafter MSFR) have shown that S(i?) decreases approximately as 
with j3 = 1 ± 0.2 to R < l/i _1 Mpc, which somewhat flattens to /3 = 0.6 ± 0.4 beyond 



l/i^Mpc to lO/i^Mpc using galaxy samples of the Sloan Digital Sky Survey (jYork et al 



1 Eq. (5) is the Taylor expansion to the first order in distortion. Menard et al. (2003) examined the accuracy of 
the formula with the conclusion that the higher order terms contribute by about 10-15 percent to magnification. 



-3 - 



2000; lAbazaiian et al.l 120051 ) . They used the correlation of magnification signal of background 



quasars with foreground galaxies. They also showed that the surface densit y profile derived from 



the magnification agrees with that from the tangential shear of galaxies (j Sheldon et al.l 120041 ) . 
Mandelbaum et al. (2006) have also derived the surface mass density profile for the luminous red 
galaxy (LRG) sample, which lies parallel to that of MSFR for R < l/i~ 1 Mpc, however, with their 
amplitudes about 2.5 times larger than that of MFSR. The difference of the surface mass density 
of the two profiles diminishes as the distance and the two agree at R > l/i _1 Mpc, albeit with 
substantial errors in the measurements. 

In this paper we study the distribution of matter around galaxies with the aid of an extensive 
use of TV-body simulations for the CDM dominated universe without baryons, which have been 
convergent among different simulations to a sufficient accuracy. We employ iV-body simulations to 
unfold contributions from multiple galaxies and to interpret the finding from gravitational lensing, 
but also to extrapolate the physics to the region beyond the reach of the observation after verifying 
that the iV-body simulation indeed describes the observed surface mass density profile with an ac- 
curacy sufficient to us. We note that Hayashi & White (2008) also used large numerical simulations 
to study the dependence of the galaxy /halo-matter cross-correlation amplitude on the galaxy mass 
and the halo mass. We expect that the state-of-the-art iV-body simulation for dark matter gives 
information that is sufficiently reliable for the quantity where the observations would not directly 
give the information, if the simulation is appropriately constrained by the observation. This gives 
insight as to the distribution of the matter beyond the observation, and would tell us concerning 
the more global distribution of matter in the Universe considerably away from galaxies. 

In section 2 we describe the iV-body simulation. In section 3 we present the surface mass 
density profile in reference to the observation. We discuss the mass distribution beyond the virial 
radius and notify that a caution is necessary as to the aperture when the mass of galaxies is referred 
to in section 4. Section 5 is a summary of our analysis. 



2. Cosmological iV-body simulations 



We use the parallelised iV-body simulation code Gadget-2 (ISpringel et al.ll200ll : ISpringelll2005l ) 
in its Tree-PM mode. Baryons are not included. We assume the flat universe with the cosmological 
parameters O m = 0.258 and h = 0.719, the Hubble constant in units of 100 km s -1 Mpc^ 1 ; for 
other parame ters n, 9 = 0.963 and <7s = 0.796, following the parameters derived in the WMAP 



5-year result (jKomatsu et al.l 120091 ) . The gravitational softening parameter is chosen to be e = 
lO/i^kpc. We take 1024 3 particles in a box of comoving 200/i 1 Mpc on a side. The mass of a 
dark matter particle is 5.34 x 10 8 }i~ 1 Mq, so that we are able to identify halos with mass a few 
times 10 10 /i _1 Mq. We set the initial redshift at zi = 50 and generate the initial condition using the 
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second-order Lagrangian perturbation theor^fS 



)Coccimarro 



19981 ). which has the advantage that it 



is more accurate to generate initial conditions so that one can s et them at an e p och la ter than with 
the conventional Zeldovich approximation; see ICrocce et al.l (|2006l ); I Jenkins! (120101) . The initia l 
matt er power spectrum at z = z% is computed using the CAMB code ([Lewis. Challinor fc Lasenbv 
20001 ). 



Haloes are identified in two-steps. We select candidate groups of dark matter particles using 
the friend-of-friends algorism, where we take the linking length to be b = 0.2. We then apply the 
spherical overdensity algorism to the candidate groups that are identified by friend-of-friends. We 
choose groups that contain at least 100 particles. We follow the conventional operational definition 
for the pseudovirial radius (r v ): when the particles enclosed within some radius give the average 
mass density 200 x p cr jt (z) , they form a halo, the mass of which is given by the sum of those particles 
and is referred to as the virial mass M v . Taking the minimum number of particles in a halo to be 
100, i.e., the minimum mass of the halo is 5.3 x 10 l0 h~ l M & , we identify 229,804 haloes. Table 1 
shows the number of haloes we identified in the simulation at z = 0.36, which is the average redshift 
of the galaxy sample used by MSFR. 

We confirmed that the halo mass function agrees with modern simulations at a good accuracy, 
e.g., with that by Crocce et al. (2010), but also with the analytic formula of Sheth & Tormen 
(1999) with a known slight deviation as reported by the above reference. We fit all groups with 
the NFW profile and compare the concentration parameter c = r v i T /r s (r s is the scale length of 
the NFW profile) as a f unction of halo mas s with ot her modern simulation s. We confirmed that 
our c agrees with that of Maccio et al. (j2008l ) (see also Bullock et al. ( 2001 )) within a few percent 
level. The simulations are well converged and the difference in simulations and analysis algorisms 
is tolerably small. 

The scale range that concerns in this paper is 20/i _1 kpc or larger, where the contribution from 
baryons modifies little the surface mass density profile. Hence, we confine ourselves to the CDM 
universe without baryons. 



2 The parallelised second-order Lagrangian perturbation theory code ( Nishimichi et al. 20091 ) is provided by T. 
Nishimichi. 



Table 1: Haloes identified in our simulation at z = 0.36. 

mass [h~ l Mo\ number number fraction [%] mass fraction [%] 

5 x 10 10 - 5 x 10 11 189,578 82^5 5^9 

5 x 10 11 - 5 x 10 12 36,402 15.8 8.25 

5 x 10 12 - 5 x 10 13 3,637 1.6 7.64 

5 x 10 13 - 5 x 10 14 187 0.08 3.13 
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3. Surface mass density profiles around haloes 

We calculate the mean surface mass density around haloes in the following way We shift 
the entire box so that it is centred on a halo being considered. We consider a beam and project 
the mass of simulation particles within comoving 100/i _1 Mpc widths to the beam. This gives the 
surface mass density around the beam at a specific projected distance. We then take the average 
of the beams around all haloes to obtain the projected surface mass density profile. 

Figure [T] compares our simulation for the surface mass density at z = 0.36 with the observa- 
tional result of MSFR which was derived at the same mean redshift. The abscissa is the physical 
distance at this redshift from the centre of the halo. The simulation is represented by a bunch of 
thin curves for 100 haloes randomly chosen from the 9970 haloes that have a virial mass larger 
than 2 x 10 12 /i~ 1 Mq (this choice is discussed below). The maximum halo mass in our simulation 
is 5 x 1O 14 /i _1 M0. The thick solid curve is the average over the entire sample above the mass 
threshold. The observational data (ticks with error bars) are taken from MSFR. The simulation 
agrees with the observation very well, in both overall shape and amplitude up to lO/i^ 1 Mpc. 

The dashed curve shows the contribution from the central haloes which are truncated at 
their virial radii, the so-called one-halo term. The figure shows that the mean surface density at 
R < 200/i~ 1 kpc is dominated by the one-halo term. We note that our column integrates over all 
particles along the line of sight and it would receive the contribution not only from the tail of 
neighbouring haloes but also from possible 'unbound' particles away from galaxies. 

To separate the contributions of bound and unbound particles, we remove the particles beyond 
the virial radii of all haloes. Then we recompute the surface mass density in the same way as 
described above, but including haloes with mass below our virial mass threshold set above. This 
is shown in Figure [2] with the dash-dotted curve, which stands for the contributions from bound 
particles in all haloes. The data, the thick solid and the dashed curves are the same as in Figured) 
The figure shows a substantial difference between the solid and dash-dotted curves, meaning that 
particles bound in haloes contribute only by 1/3 the total surface mass density beyond the distance 
of ~ 500/i _1 kpc from the galaxy. We remember that the total surface mass density derived in the 
simulation agrees with that estimated from gravitational lensing. Roughly 2/3 the surface mass 
density at such a distance is due to particles residing beyond the virial radius of any galaxies, i.e., 
gravitationally 'unbound' particles. 

In our argument we set the cutoff on the lowest mass of haloes, somewhat arbitrarily, at 
2 x 10 12 /i _1 Mq. This choice should of course affect the surface mass density profile. We show in 
Figure [3] the mass density profiles with different threshold mass: 5 x 10 n ,2 x 10 12 ,5 x 10 12 ,2 x 
10 13 ,5x IO^^Mq. We see that the surface mass density increases as the lower cutoff mass (virial 
mass) M\ ow increases for small radii, say, R < lft,~ 1 Mpc where the one halo term is the significant 
contributor. The result of the simulation shows S(i?) oc MJ for the distance scale lOkpc to a few 
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hundred kpc: it is summarised, e.g., at 120ft, 1 kpc as 

£ ~ 100(M low /8 x lO 12 /i~ 1 M ) 2 / 3 /iM Q (pc)- 2 



(6) 



With the singular isothermal sphere we have the surface mass density at the projected distance R, 



(7) 



and a similar relation holds with the NFW profile albeit in the limited distance range. With our 
mass function the average mass (M v ) ~ 5Mi ow , so that our fitting formula eq.([6]) gives for a given 
S(i?) a halo mass 2.5 times smaller than the model with the singular isothermal sphere. 

This sample threshold dependence exp lains the difference in the surface mass density distri- 
butions between iMandelbaum et al.l (|2006l ) . which are also plotted in Figure EJ and MSFR. The 
former gives Ti(R) larger by a factor of 2.5 than MSFR within a few hundred kpc range. Equation 
(6) then indicates that the threshold mass of the LRG sample of Mandelbaum et al. (2006) is ap- 
proximately by 4 times more massive than that of MSFR, who used the main galaxy sample. The 
difference seen in among different halo masses diminishes for a large R, where the one halo 

term no longer dominates but E(-R) is contributed by ne ighbouring haloes and unbound particles, 
as seen by comparing S(i?) of IMandelbaum et al.l (B) and MSFR. @ (Neighbouring haloes are 
more likely to be those of normal galaxies rather than LRGs.) 

Figure [3] indicates that the mean surface mass density profile around galaxies measured by 
MSFR is reproduced well when the threshold mass M ]nwr is set to 2 x W^h^M ^, which is approx- 
imately the mass of the L* galaxy 1.2 x 1O 12 /i -1 M (jFukugita fc Peebles! 12009 ): hence our default 
cutoff is chosen. MSFR used the galaxy sample 17 < i < 21. If the effective cutoff of th e sample 
would be around i ~ 20, sample's threshold is around the L* galax\F( Blanton et al. 2001 ). 



It is less ambiguous to estimate the threshold mass from T,(R) itself. With the threshold 
2 x 10 h~ l MQ, of the simulation is in good agreement with the data. If we would take the 

threshold 5 x 10 11 /i _1 M , which roughly corresponds to the i ~ 21 mag threshold, of the 

simulation would lie somewhat too low. This consideration suggests that the surface mass density 
can be used to characterise the mass of the galaxy sample, which otherwise is difficult to estimate 
accurately. 

In the limited range of the distance scale we are studying, 15kpc < R < 200kpc both NFW 
and singular isothermal sphere give mass profiles that are very similar and both give good fits to 
the data, so that we are not able to distinguish between the two profiles in this range. 



3 The LRG data show a slightly larger amplitude at R ~ lh 1 Mpc, reflecting gre ater bias. The ove r all fea tures 
of the surface density profiles for different galaxy/halo masses are studied in detail by Havashi fc White] (2008). 

4 At a more accurate level we must consider that L* luminosity corresponds to 10 10 ' 6 Lq in the i band, while it 
does to lO lo,5 I/0 in the r band. L* luminosity is not physically well defined. We also note that there is a significant 
uncertainty in the mass-luminosity relation associated with morphology of galaxies. 
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4. Mass distribution beyond the virial radius 



We see above that the substantial amount of matter in the Universe resides outside the pseu- 
dovirial radius of galaxies. In order to examine whether the distribution of unbound particles 
is organised, we calculate the total amount of mass encircled with the pseudovirial radius, and 
then expand the encircling radius by a factor of a: we denote it in units of the critical density 
as f^haio extended («)• The pseudovirial radius is defined for each galaxy as 200 times p cr it- For a 
large a two haloes start overlapping, and in such cases we count the amount of material only once, 
avoiding double counting. Figure [H shows the fraction Ohaio extended (aO/fi m as a function of a, 
which is described well with 0.23 In a + 0.22 consistent with the NFW profile in its p ~ r~ 3 regime. 

If the mass distribution of a galaxy were spatially bounded, it would not matter for the mass 
estimate of galaxies whatever the encircling radius is taken in so far as it is large enough. We find, 
however, this is not the case. The encircling radius must be carefully specified, for instance, to 
calculate the mass to light ratio to be quantitatively meaningful. 

The figure shows that the matter within the pseudovirial radius (a = 1) is only 0.25 times the 
total matter if the sample of galaxies is set above the mass thereshold 10 11 h~ 1 MQ, above which 
the galaxy mass (and hence M/L) is observationally estimated; see below. This fraction (and the 
curve in the figure) depends on the threshold mass of sample galaxies, as ~ — log M\ ow in so far 
as M\ ow is taken below the break mass of the mass function. If the threshold is taken at a larger 
mass, this fraction decreases faster, mainly because the number of galaxies included in the sample 
decreases. Note that the sample threshold must be taken small enough to estimate the mass density 
of the universe so that the contribution of sub-threshold galaxies to the global quantity (in practise 
luminosity density) is not substantial. We confirm that this mass fraction agrees with the mass 
obtained by integrating over the mass function. The figure also shows that this fraction becomes 
close to 0.5 if 3 times the pseudovirial radius is taken to encircle clustering. Almost entire matter 
(> 90%) is included only with a > 20. 



Now it may be appropriate to examine the case discussed by iMcKay et al.l (120011 ). who es- 
timated the mass of galaxies from gravitational lensing shear, as we quoted earlier. They used 
the spectrosco pic sample of SPS S, which means the limiting magnitude accurately defined at 

—20.6 for the median redshift of 



r iimit = 17.8 (IStrauss et al 
the sample z 



2002). It corresponds to M r 



0.10 with the median K correc tion, or to 0.41L*. If we invoke an empirical scal- 
ing relation M oc I? (jMandelbaum et al.ll2006l ) as consistent with the Faber- Jackson relation, we 
estimate the limiting mas s of the sample M\ ow sa 2 x 10 n /i _1 Mq, taking Mmw ~ 2.3 x 1O 12 M 
(jFukugita &: Peebles! 120061 ) and Lmw = 1-4L* for the Milky Way to normalise the M — L relation. 



This can also be verified from the gravitational lensing measurement itself. IMcKay et al.l (120011 ) 
obtained E(R) w 2.5(i?/l/i~ 1 Mpc)~ a8±a2 ft,M pc~ 2 . Consulting with our £(i?) relation with the 
sample threshold mass as a parameter (see Figure [3|), we estimate M\ ow ~ 1.5 x 10 11 /i^ 1 Mq in 
agreement with the photometric estimate. 
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With this threshold we estimate the mean virial mass of the spectroscopic sample approxi- 
mately M v ^ 1.5 x 10 12 /i _1 M Q using the mass function obtained in our simulation. This value is 



compared with the estimate M v ss 2.6 ± 1.3 x 10 12 h 1 M of iMcKav et al.1 (|200ll ) for their sam 



pie, where the error includes the variance associated with the galaxy morphology. This shows the 



consistency of the two estimates, those by lMcKay et al.1 ()200ll ) and ours. 



Our simulation gives the average pseudovirial radius over the sample with the lower mass cutoff 
Mi ow as r v = lOOft, _1 kpc(Mi ow /lO 11 /i _1 M0) ' 29 , so that our estimate of the virial radius TZ for the 



sample with M iow m 1.5 x 10 11 h~ 1 M Q is r v = 120/i _1 kpc. IMcKav et al.l (|200ll ) claim that they 
measured the M/L within 260ft, -1 kpc, which is a ~ 2.2 times the mean virial radius of the sample. 
Figure d] tells us that for a = 2.2 approximately 0.4 times the total Vl m should be included in their 
estimate. This leads us to the global matter density from galaxies to be 

n m = (0.13 ±0.05)/0.4 « 0.32, (8) 

which is consistent with the global value of Eq.©. Although our estimate presented here is ad- 
mittedly crude, this agreement indicates that the mass beyond the pseudovirial radius we inferred 
here is probably broadly correct. We should underline the importance of the estimate of the radius 
(e.g., with respect to the virial radius) when the mass or the mass to light ratio is presented for 
galaxies, since galaxies are extended objects without definite boundaries. 

We next calculate the volume occupied with extended haloes as a function of a. In Figure [5] 
we plot the fraction of mass contained in the extension of haloes Sl^aio extension ( a ) /Vl m as a function 
of the fraction of the volume occupancy Vhaio extension (oOAtotab where Hotai is the total simulation 
volume. The numbers beside the symbol indicates the multiplier a. Vh a i extension is approximately 
proportional to (aR) 3 up to overlaps of haloes at a large a. This figure shows that the plot is given 
nearly by a straight line, indicating that 

^halo extension {a) /Vt m ~ [Vhalo extension («) /Vtotal] ' (9) 

Let us recall if the mass distribution is random throughout the entire volume we expect fihaio extension cx 
^haio extension- The figure shows that there is no symptom of a conspicuous break of the curve at 
any a, meaning that the distribution of unbound matter is all organised to a distance significantly 
away from the galaxy, without leaving a significant amount of material in the intergalactic space. 
This power means that the matter density behaves as 

p m ~ r~ 2A , (10) 

beyond the virial radius, as long as halfway to the neighbouring galaxy. This drops faster than 
the isothermal profile but is numerically consistent with the tail of the NFW profile in the range 
r/r s = 5 — 100 with r s the NFW scale radius. Remembering that r v /r s ~ 5 for galaxies and the 
typical spacing between the galaxies is tq ~ 5/i -1 Mpc, the relevant range matches the inter-galaxy 
distance. 
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Since galaxies have no clear edges, it is not appropriate to define their 'total' mass. If we 
include the mass in the tail of the galaxy halfway to the neighbouring galaxy, the mass of galaxy 
increases approximately by a factor of 2 so that the effective M/L becomes ~ 35O/IM0L" 1 . When 
multiplied by the luminosity density L = 2.2 x 10 8 /i.LqMpc _3 , we would obtain 

p m = L x (M/L) « 0.3/9 cri t, ft m «0.3. (11) 

We remark that clusters serve as a good natural integrator of the mass of galaxies. The mass of 
clusters is usually estimated at the radius, say, rsoo which is far beyond the virial radii of individual 
member galaxies. Therefore mass of galaxies resided in the tail is largely integrated in the estimate 
of the cluster mass. This explains the reason why M/L of clusters reaches 300—400, significantly 
larger than the estimates for individual galaxies: the mass here includes the mass present in the 
tail of galaxies. This explains the reason why we are arrived at the correct global mass density 
of the Universe, if the cluster M/L is used to multiply on the luminosity density of field galaxies 
instead of M/L of individual galaxies, though this is apparently an incongruous treatment. 



5. Summary 

We showed that the state of the art iV-body simulation of dark matter based on the ACDM 
model gives an excellent description of the surface density profile of the mass distribution around 
galaxies, which has been explored up to the lO/i" 1 Mpc scale from the galaxy mass correlation 
function using weak gravitational lensing analysis applied to large modern galaxy samples. The 
surface profile thus derived is consistent with r -1=t0 ' 2 up to l/i _1 Mpc and somewhat flattens to 
r -o.6±o.4 bgycmol this radius. The latter power is consistent with that of the two point correlation 
function of galaxies. The galaxy mass correlation function measures the mass distribution within 
the galaxy halo in a short distance scale (smaller than a few x 100 kpc scale) and reflects the 
galaxy distribution beyond this radius. The two distributions are similar and match each other 
with a slight break. 

The amplitude of the surface mass density profile depends on the galaxy sample, in particular 
on the lower mass cutoff applied to the sample. Hence the amplitude serves us to infer the properties 
of the galaxy sample in a self-consistent way. 

We showed that the galaxy has no clear edges in the dark matter distribution, unlike luminous 
matter, which should be bounded by the cooling radius. The distribution is extended beyond 
the virial radius in an organised way halfway to the neighbouring galaxy, so that the Universe 
is filled with the material associated with tails of galaxies, and we then call the peaks of the 
matter distribution galaxies. Inter-galactic space is filled with matter. Tails of galaxies extend to 
great distances without cutoff, whereas luminous component of galaxies have definite cutoff radius 
corresponding to the cooling radius. 

About half the matter in the Universe is gravitationally unbound at z ~ 0. Its distribution, 
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however, is never random or uniform, but is well organised in a way to be consistent with the 
tail of galaxies with the mass density roughly p ~ r~ 2A . Half the matter is present in the tail of 
galaxies beyond the pseudovirial radius. This explains the gap in the estimate of mass density of 
the Universe between the global value and the value obtained by adding the contributions from 
matter bound to individual galaxies: with this extende d matter distribution the m atter entry closes 
in the mass inventory, which has been left unclosed in iFukugita fc Peeblesl {200M . 



The observables derived from gravitational lensing lend us to characterise the sample in terms 
of the mass with physical means. The agreement between photometrically inferred characteristics 
and physically derived ones is, if not perfect yet, nearly satisfactory. Gravitational lensing would 
be used to characterise the sample with physical methods. 
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Fig. 1. — The surface mass density profile as a function of the distance from the centre of galaxies. 
The bunch of thin solid curves represent the profiles for 100 randomly picked-up haloes from 9970 
haloes with mass > 2 x 10 12 /i~ 1 Mq genarated in the ACDM iV-body simulation. The thick solid 
curve is the mean of all haloes above the mass threshold. The dashed curve shows the contribution 
from the one-halo term. The data with error bars are the observational estimate using galaxy mass 
correlation function deduced from gravitational weak lensing for quasar brightness in MSFR. The 
abscissa is the physical distance at z = 0.36, which is the average redshift of the weak lensing 
observation. 
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Fig. 2. — The mean surface mass density profile as a function of the distance from the centre of 
galaxies. The thick solid curve is the mean of all haloes above the mass threshold. The dash-dotted 
curve represents the contribution from particles bound to haloes, i.e., particles that reside within 
the virial radius of all haloes. The data with error bars are the observational estimate by MSFR 
as in the previous figure. 
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Fig. 3. — The mean surface mass density profile as a function of the distance from the centre of 
galaxies for galaxy samples with the differ ent threshold halo mass (y irial mass). The open squares 
which represent the data for LRG given by iMandelbaum et al.l (120061 ) are added to the MSFR data 
for the SDSS main galaxy sample shown with solid diamonds with error bars. 
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Fig. 4. — Fraction of mass contained in the sphere centred on individual haloes with radius aR v i r , 
where i? v i r is the pseudovirial radius and a is the multiplier represented in the abscissa. The plot 
is for the sample with the threshold halo mass 1 x 10 11 M Q . The solid curve is 0.23 In a + 0.22 given 
in the text. 
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Fig. 5. — Mass fraction borne by extended haloes with the radius r < aR vir , Q\ VA \ extension («)/^m 
versus the volume fraction occupied by the extension of virial spheres Vhaio'extension(oO/^tot- The 
numbers shown beside the symbol is the multiplier a. The solid line indicates O n aio extension/^m ~ 

(^halo extension 



